Pseudorapidity distributions of charged particles in asymmetric collisions using Tsallis 
thermodynamics* 


J.Q. Tao,! H.B. He, H. Zheng,” Ý W.C. Zhang,” X. Liu,? L.L. Zhu, and A. Bonasera* 6 


'Key Laboratory of Quark & Lepton Physics (MOE) and Institute of Particle Physics, 
Central China Normal University, Wuhan 430079, China 
>School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710119, China 
4 Institute of Nuclear Science and Technology, Sichuan University, Chengdu 610064, China 
‘Department of Physics, Sichuan University, Chengdu 610064, China 
Cyclotron Institute, Texas A&M University, College Station, TX 77843, USA 
®Laboratori Nazionali del Sud, INFN, 95123 Catania, Italy 


The pseudo-rapidity distributions of the charged particles produced in the asymmetric collision systems p+Al, 
p+Au and 3He+Au at /8nNn = 200 GeV are evaluated in the framework of a fireball model with Tsallis ther- 
modynamics. The fireball model assumes that the experimentally measured particles are produced by fireballs 
following the Tsallis distribution and it can effectively describe the experimental data. Our results as well as 
previous results for d+Au collisions at snn = 200 GeV and p+Pb collisions at \/s~n = 5.02 TeV validate 
that the fireball model based on Tsallis thermodynamics can provide a universal framework for pseudo-rapidity 
distribution of the charged particles produced in asymmetric collision systems. We predict the centrality de- 
pendence of the total charged particle multiplicity in the p+Al, ptAu and ?He+Au collisions. Additionally, 
the dependences of the fireball model parameters (Yoa, YOA, Ca and g4) on the centrality and system size are 


studied. 
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I. INTRODUCTION 


High-energy heavy-ion collisions provide a unique way to 
understand the origin of the universe. However, their pro- 
cesses cannot be directly observed in experiment. We can 
only study the collision process indirectly by analyzing the 
properties of the final particles produced in the collisions. The 
pseudo-rapidity distribution of charged particles is one of the 
important experimental observables. The study of this ob- 
servable could lead to a better understanding of the properties 
of the particles produced in the collisions, the particle produc- 
tion mechanism and so on. There have been numerous works 
in previous studies using different models, such as HIJING 
[1], AMPT [2-4], EPOS-LHC [5], a multi-source thermal 
model [6, 7], a new revised Landau hydrodynamics model 
[8], a 1 + 1 dimensional hydrodynamics model [9, 10], a dy- 
namical initial state model coupled to (3 + 1)D viscous rel- 
ativistic hydrodynamics [11] and so on, to analyze the exist- 
ing experimental data of pseudo-rapidity distributions of the 
charged particles [12-30]. Although these models are based 
on different physical ideas, valuable physical information on 
the collision process has been extracted and learned. 

Recently, a fireball model based on Tsallis thermodynam- 
ics was utilized to analyze the pseudo-rapidity distribution of 
the charged particles measured in high-energy heavy-ion col- 
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lisions [31-33]. In our previous works [32, 33], we used the 
fireball model to study the pseudo-rapidity distributions of the 
charged particles produced in p+p(p) collisions for energies 
ranging from ,/syn = 23.6 GeV to 13 TeV and A+A colli- 
sions at the RHIC and LHC, and extended the fireball model 
to the asymmetric collision systems, i.e., d+Au collisions at 
/snn = 200 GeV and p+Pb collisions at syyn = 5.02 
TeV, by considering the asymmetric collision geometry con- 
figuration. In this paper, we utilize recently published data 
from the PHINEX Collaboration [30] at RHIC to systemat- 
ically study the pseudo-rapidity distributions of the charged 
particles produced in asymmetric collision systems, includ- 
ing p+Al, p+Au and ?He+Au collisions at V/snn = 200 
GeV. We also predict the total multiplicities of the charged 
particles from the fireball model and study their centrality de- 
pendence. Further, we analyze the centrality and system size 
dependencies of the fireball model parameters obtained from 
the pseudo-rapidity distributions of the charged particles. 


The paper is organized as follows. In section II, the fireball 
model with Tsallis thermodynamics is briefly introduced. In 
section III, the fitting results of the fireball model and the to- 
tal charged particle multiplicities extracted from the fireball 
model are shown. The dependences of the model parame- 
ters on the centrality and size of the collision systems are also 
presented. A brief conclusion is drawn in section IV. 


Il. THEORETICAL DESCRIPTIONS 


In the self-consistent Tsallis thermodynamics, the Tsal- 
lis distribution is proposed as a generalization of the Boltz- 
mann-Gibbs distribution [34]. To describe the transverse mo- 
mentum spectrum of particles, the Tsallis distribution is writ- 
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s7 where g is the particle state degeneracy, V is the volume, 
ss Mp = \/mé~ + p> is the transverse mass and mp is the parti- 
sə Cle rest mass, y is the rapidity, q is the entropic factor, which 
e Measures the non-additivity of the entropy [34, 35], u is the 
61 chemical potential and T is the temperature. The Boltzmann 
e2 distribution is recovered when q = 1. We take u = 0 because 
63 the multiplicities of 7+ and 77 are equal and they are the 
64 Majority of particles produced in the collision systems con- 
e sidered. For the middle rapidity y ~ 0, Eq. (1) can be rewrit- 
es ten as 
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es The parameters q and T are extracted from the experimental 
69 transverse momentum spectrum of the particles. 

70 In the fireball model with Tsallis thermodynamics [31- 
71 33], the particles measured in the experiment were produced 
72 by fireballs following the Tsallis distribution Eq. (1). The 
73 density distribution of these fireballs in the rapidity space is 
74 v(yf), Where yy is the rapidity of the fireball. Therefore the 
75 transverse momentum spectrum of particles can be written as 
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J 73 where N is the total particle multiplicity and A is the normal- 
79 ization constant such that 
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81 Sometimes, the experimental data are measured in the 
s2 pseudo-rapidity 7 space. To describe the experimental data 
es “> we substitute the relation between rapidity and pseudo- 
s4 rapidity [36] 
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ss into Eq. (3) and integrate the transverse momentum in the 
s7 equation to obtain [32, 33] 
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Because of the term of 4/1 — eee Eq. (6) cannot be 


analytically integrated over pr and it is done numerically. 
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Fig. 1. (Color online) The pseudo-rapidity distributions of the 


charged particles produced in p+Al collisions at ,/s~nv = 200 GeV 
for different centralities. The symbols are experimental data taken 
from [30]. The curves are the results from Eqs. (6) and (8). 
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Fig. 2. (Color online) Same as Fig. 1, but for p+Au collisions at 
y SNN = 200 GeV. 
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Fig. 3. (Color online) Same as Fig. 1, but for 7He+Au collisions at 
/SNN = 200 GeV. 


asymmetric g-Gaussian functions [33], 
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where Yoa(A) and Caça) are the centroid position and width 
of fireball distribution in the direction of the light (heavy) nu- 
cleus beam, respectively. The normalization of Eq. (8) is 
handled by the normalization constant A in Eq. (3). x is the 
parameter to characterize the extent of asymmetry which was 
first proposed in our previous work [33]. In this work, we 
take q’ = q same as in [31-33]. A representative figure of Eq. 
(8) is shown in the Appendix A with the parameters obtained 
for the p+Al collisions at 0-5% centrality and \/snnw = 200 
GeV. 
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Ill. RESULTS AND DISCUSSION 


Because the data of the transverse momentum spectra of 
the charged particles produced in p+Al, p+Au and 3°He+Au 
collisions at ,/sy7;7. = 200 GeV have not been released yet, 
the data of the transverse momentum spectra of 7° produced 
from these collisions are obtained from [37] in this study. Us- 
ing Eq. (2) by taking g = 1 and V as a free parameter, the 
parameters T and q for the Tsallis distribution are extracted 
and listed in Table 2 in the Appendix B. A representative fig- 
ure of transverse momentum spectra of 7° for p+Al collisions 
at \/Svn 200GeV is shown in Appendix C. It is worth 
noting that the transverse momentum spectrum of 7° is very 
similar to that of 7* at ,/syyN = 200 GeV for a given col- 
lision centrality, the temperature parameter T extracted from 
n° should reasonably characterize the property of the colli- 
sion system. We take the parameters T and q from the closest 
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Fig. 4. Total charged particle multiplicities produced in the p+Al, 
p+Au and 3He+Au collisions at V/snw = 200 GeV versus the col- 
lision centrality c. The squares are the fireball model results. The 
lines are the fitting results. The fitting function is from [14] and 
specified in the legend. 


centrality when the centrality of the particle transverse mo- 
mentum spectrum and the centrality of the charged particle 
pseudorapidity distribution are not the same. These two pa- 
rameters and the fireball model with Tsallis thermodynam- 
ics, Eqs. (6) and (8), are then utilized to study the pseudo- 
rapidity distribution of the charged particles produced in the 
collisions. The corresponding values of x in Eq. (8) are also 
listed in Table 3 in the Appendix B. 

In Figs. 1, 2 and 3 the results of the pseudo-rapidity dis- 
tributions of the charged particles from the fireball model 
with Tsallis thermodynamics for different centrality bins in 
pt+Al, p+Au and ?He+Au collisions at /s8nn = 200 GeV 
are shown. The fireball model effectively describes the exper- 
imental data within the errors. Notably, the data quality of the 
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Fig. 5. (Color online) The centrality dependence of model parameters Yoa, YoA, Oa and o4 in p+Al, p+Au and 2He+Au collision at \/snn = 
200 GeV. The lines are the linear fit results to guide the eyes. 
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Fig. 6. (Color online) The collision system size dependence of model parameters you, yoA, Ca, 7A for p+p, ptAl (0-5%), ptAu (0-5%), d+Au 
(0-20%), ?>He+Au (0-5%) and Au+Au (0-6%) collisions at ,/sn jn = 200 GeV. The parameters for the p+p, d+Au and Au+Au collisions are 
taken from Ref. [33]. 


137 pseudo-rapidity distributions of the charged particles is notas 142 and affects our analyses of the fireball model parameters ver- 
138 good as that for d+Au collisions at \/s;y~7 = 200 GeV shown 143 sus collision centrality and collision system size to some ex- 
139 in [33], i.e., in terms of larger errors, a fewer number of data 144 tent later in the following. The pseudo-rapidity distribution 
140 points as well as a lower pseudo-rapidity coverage, which 145 of the charged particles for centrality 5-10% is lower than the 
141 leads to larger uncertainties to the fireball model parameters 14 case for centrality 10-20% in some pseudo-rapidity regions 


147 for 3He+Au collisions which is observed in Fig. 3. A larger 
x at centrality 5-10% compared with the others for 7He+Au 
collisions is also observed in Table 3. We emphasize that the 
same fitting protocol is applied for all the pseudo-rapidity dis- 
tribution data of the charged particles. Because the collision 
system is asymmetric, the pseudo-rapidity distribution of the 
charged particles has significant forward/backward asymme- 
try. Fewer particles are produced in the direction of the light 
nucleus (p, >He) beam compared to the heavy nucleus (AI, 
Au) beam. As the d+Au collision system at ,/svn = 200 
GeV and the p+Pb collision system at \/svn = 5.02 TeV 
we studied in [33], the pseudo-rapidity distributions of the 
charged particles produced by these collision systems also 
become more symmetric from the central to peripheral col- 
lisions. It is because the peripheral collisions for asymmetric 
collision systems are more similar to the symmetric p+p col- 
lisions according to collision geometry. 
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Table 1. Results of the linear fits are shown in Fig. 5. The c repre- 
sents the centrality. 


System (GeV) Yoa Oa Yoa OA 


p+Al 200 
p+Au 200 


2.32c+1.33 -0.19c+1.57 -2.13c+2.04 1.51c+1.32 
0.61c+1.16 0.56c+1.09 -0.69c+2.19 -0.27c+1.95 


3Het+Au 200 -0.60c+1.96 2.03c+0.26 0.32c+1.63 -1.21c+2.32 


We then evaluate the centrality dependence of the total 
multiplicities of the charged particles produced in these colli- 
sion systems. Integrating Eq. (6) over the 77 € [—10, 10] we 
obtain the total multiplicity of the charged particles for each 
centrality from the fireball model. Because the corresponding 
experimental data are not yet available, we only analyze the 
results extracted from the fireball model and treat them as pre- 
dictions. Figure 4 shows the total multiplicities of the charged 
particles calculated from the fireball model versus collision 
centrality c. c = 0 represents the most central collisions and 
c = 1 represents the most peripheral collisions. It can be ob- 
served that the fitting function taken from [14] can effectively 
describe the centrality dependence of the total multiplicities 
of the charged particles. As the centrality changes from the 
178 central to peripheral collisions, fewer charged particles are 
179 produced. 

We also analyze both the centrality and system size depen- 
181 dence of the parameters (Yoa, YoA, Ca and o4) of the fireball 
182 model. In Fig. 5 the dependence of the fireball model pa- 
183 rameters on the collision centrality in the p+Al, p+Au and 
isa °>He+Au collisions at /8NN = 200 GeV are shown. Inspired 
185 by the linear relation of the centrality dependence of the fire- 
186 ball model parameters for d+Au collisions at ,/svn = 200 
187 GeV and p+Pb collisions at ,/syn = 5.02 TeV shown in 
188 Fig. 12 of our previous work [33], the linear fittings are per- 
189 formed to guide the eyes in Fig. 5 and the fitting functions are 
190 listed in Tab. 1. The negative and positive slopes of the linear 
191 fittings for the fireball model parameters (yoa, Yoa and ca) 
192 versus centrality are similar to those of the d+Au and p+Pb 
193 Collisions in [33]. In the following discussion, the results for 
194 d+Au and p+Pb are obtained from [33]. It can be observed 
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195 that the slopes of the linear fittings for parameter yoa ver- 
196 SUS Centrality are positive and the corresponding slopes of the 
197 linear fittings for parameter yoa versus centrality are nega- 
198 tive for p+Al, p+Au and d+Au collisions at /syn = 200 
199 GeV. However, the slopes of the linear fittings for parame- 
200 ters (Yoa, Yoa) reverse the signs respectively for 3He+Au col- 
21 lision at ,/s~nw = 200 GeV compared to the above men- 
22 tioned cases. The slopes of the linear fittings for parameters 
3 (Yoa, Yoa) are positive for p+Pb at ,/svn = 5.02 TeV. It can 
also be observed that there is a universal trend with increas- 
ing centrality for parameter 0,4) except the lightest colli- 
sion system p+Al, i.e., 7, increases with increasing centrality 
in the direction of the light nucleus beam and a, decreases 
with increasing centrality in the direction of the heavy nucleus 
beam. In the p+Al collision system, ca and ø 4 have opposite 
trends with increasing centrality to their counterparts in the 
other collision systems. These different patterns indicate the 
complex dynamics in the asymmetric collisions relevant to 
the combinations of the projectile and target as well as the 
collision energy which needs more investigations. 

Figure 6 shows the collision system size dependence of the 
fireball model parameters at \/s Nw = 200 GeV. For collision 
systems other than p+p, the parameters of the most central 
collisions are considered. In the p+p and Au+Au collisions, 
the parameters you = YoA = Yo. Ta = TA = O, Where yo and 
o are the rapidity centroid and width of fireball distribution 
in the symmetric collision system, as detailed in [33]. It can 
be deduced that when the light nucleus is p, Yoa decreases as 
the size of the heavy nucleus increases, whereas yo 4 shows 
the opposite trend. This indicates that a larger heavy nucleus 
has stronger stopping power for p. When the heavy nucleus 
is Au, Yoa increases as the size of the light nucleus increases, 
whereas yoa Shows the opposite trend. This means that a 
larger light nucleus is more difficult to be stopped by Au but 
has a stronger stopping power for Au. For parameters o, 
and g4, no conclusive patterns are observed. We expect that 
231 More discussions can be added when the data quality of the 
232 pseudo-rapidity distributions of the charged particles is im- 
233 proved by experimentalists. These phenomena manifest the 
complex dynamics in the asymmetric collisions. 
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IV. SUMMARY 


235 


In this paper, we studied the pseudo-rapidity distributions 
of the charged particles produced in p + Al, p + Au and 
3He + Au collisions at /snn = 200 GeV using the fireball 
model with Tsallis thermodynamics. The model can well fit 
the experimental data from asymmetric collisions. We also 
extracted the total multiplicities of the charged particles from 
the fireball model as predictions and analyzed their depen- 
243 dence on collision centrality. Notably, the data quality of the 
244 pseudo-rapidity distributions of charged particles produced 
245 inp + Al, p + Au and 3He + Au collisions at Vsnn = 200 
246 GeV affected our results to some extent. Combining our 
247 previous results of d+Au collisions at ysy = 200 GeV 
z8 and p+Pb collision at \/svn = 5.02 TeV, we analyzed 
249 the centrality and system size dependence of the fireball 
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model parameters (Yoa, YoA> Ca and o 4). Interesting patterns 
were revealed, which indicated the complex dynamics 
in the asymmetric collisions. Our results confirmed the 
conclusion made previously in [33] that the fireball model 
with Tsallis thermodynamics as a universal framework could 
also describe the pseudo-rapidity distribution of the charged 
particles produced in asymmetric collision systems. 
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Appendix A: Fireball distribution of Eq. (8) 


A representative figure of the fireball distribution of Eq. (8) 


272 using the parameters obtained for the p+Al collisions at 0-5% 
273 centrality and ,/sj nv = 200 GeV is shown. 
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Fig. 7. (Color online) Fireball distribution with the parameters ob- 
tained from p+Al collisions for 0-5% centrality at \/snn = 200 
GeV is shown. (a) The value of x is varied; (b) The value of q’ is 
varied. 


Appendix B: Parameters q, T and x 


The parameters of q and T extracted by fitting the trans- 
verse momentum spectrum of particles [37] with Tsallis dis- 
tribution Eq. (2) as well as the parameter x in Eq. (8) are 
listed. 


Table 2. Parameters q and T for the p+Al, p+Au and ?He+Au colli- Table 3. Parameter x for the p+Al, p+Au and *He+Au collisions at 
sions at snn = 200 GeV for different centralities. snn = 200 GeV for different centralities. 
System Centrality q T(GeV)  y?/NDF System Centrality x 
pt+Al 0-5%  1.092+0.004 0.127+0.012 0.129 pt+Al 0-5% 1.7600.796 
0-20%  1.098+0.004 0.115+0.012 0.151 5-10% 1.313£0.620 
20-40% 1.101+0.004 0.107+0.010 0.288 10-20% 1.247+0.596 
40-60%  1.101+0.004 0.1080.012 0.223 20-40% 1.430-+£0.840 
60-72% 1.103+0.005 0.098+0.013 0.232 40-74% 4.767+4.522 
p+Au 0-5% 1.0890.004 0.138+0.012 0.261 p+Au 0-5% 4.098+1.540 
0-20%  1.090+0.004 0.138+0.012 0.117 5-10% 2.318+0.190 
20-40% 1.095+0.004 0.127+0.012 0.122 10-20% 4.690+1.637 
40-60%  1.096+0.004 0.123+0.012 0.204 20-40% 3.173+1.196 
60-84% 1.105+0.004 0.103+0.010 0.164 40-60% 1.470+0.194 
3He+Au 0-5%  1.092+0.004 0.133+0.010 0.082 60-84% 1.092+5.980 
0-20%  1.095+0.004 0.127+0.011 0.066 3He+Au 0-5% 5.912+1.926 
20-40%  1.098+0.004 0.118+0.011 0.046 5-10% 13.055+7.709 
40-60%  1.101+0.004 0.112+0.011 0.164 10-20% 5.470+1.900 
60-84% 1.103+0.004 0.104+0.011 0.157 20-40% 5.414+1.919 
40-60% 2.210+0.957 
60-88% 1.184-£0.550 
281 Appendix C: The particle spectra fit 


22 Using Eq. (2), the fitting figure of the transverse mo- 
23 mentum spectra of 7° produced in the p+Al collisions at 
234 4/SNN = 200 GeV is shown. Similar results are obtained 
25 for the p+Au and °He+Au collision systems at /snn = 200 
286 GeV. 
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Fig. 8. (Color online) Transverse momentum spectra of 1° produced 
in p+Al collisions at ,/s~nw = 200 GeV [37]. The curves are the 
corresponding fittings using Eq. (2). 
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